Holon edge states in organic conductors (TMTTF)2X in the charge-gap regime 
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Using the bosonization technique, which separates charge and spin degrees of freedom, we study a 
possibility of formation of the holon edge states in a one-dimensional electron system with an energy 
gap in the charge sector. The results are applied to the quasi-one-dimensional organic conductors 
(TMTTF)2X. The different roles of the bond and site dimerizations in this material are discussed. 
We predict that the holon edge states should appear below the temperature of the recently discovered 
ferroelectric transition, where the inorganic anions X displace to asymmetric positions. 



When the bulk electron states have an energy gap, ad- 
ditional states may exist on the surface with the energies 
inside the gap. Such surface states were studied for a long 
time in the cases where the gap has a band-structure ori- 
gin . More recently, edge states attracted much interest 
in the cases where the gap originates from a phase transi- 
tion in the electron system. Examples are midgap states 
in the high- Tc d-wave cuprate superconductors |^ and 
chiral states in the p-wave superconductor Sr2Ru04 (see, 
for example, Ref. |^ and references therein). Edge states 
in the quasi-one-dimensional (QID) organic conductors 
(TMTSF)2X H were studied for the triplet supercon- 
ducting state H and for the magnetic-field-induced spin- 
density- wave (FISDW) state |^ , which exhibits the quan- 
tum Hall effect. In all of these case, noninteracting elec- 
trons were subject to a mean-field potential, so the prob- 
lem reduced to a single-particle Schrodinger-like equa- 
tion with a boundary. However, it is known that a one- 
dimensional (ID) system of interacting electrons can be 
described in terms of separate charge and spin sectors 
(see, for example, review Q]). An edge state in such a sys- 
tem with a gap in the spin channel was studied in Ref. fel . 
In this Letter, we use the technique developed in Ref. |g] 
(see also Ref. j9| ) to investigate another physical situation 
where a gap opens in the charge sector. The charge gap 
appears when the electron system is commensurate with 
the crystal lattice, and the umklapp interaction becomes 
relevant. This is the case for the (TMTTF)2X family of 
QID organic conductors [0, which are the best known 
examples of spin-charge separation. In these materials, 
electric transport is thermally-activated due to the charge 
gap, whereas spin susceptibility only weakly depends on 
temperature, which indicates gapless spin sector. 

The (TMTTF)2X crystals consist of parallel ID molec- 
ular chains quarter-filled with holes. The umklapp term 
appears due to dimerization of the chains. Gener- 
ally, there are two sources of dimerization, one due to 
nonequivalence of hopping integrals along the chain, and 
another due to nonequivalence of the site potentials. The 
unit cell consists of two organic molecules TMTTF and 
one inorganic anion X. The anions X are located out- 
side of the TMTTF chain, between the pairs of TMTTF 
molecules, as shown in Fig. |l|. Because of the anions, the 
bonds between the TMTTF molecules alternate between 



short and long. At high temperatures, the anions occupy 
symmetric positions, so the TMTTF molecules them- 
selves are equivalent. However, it was recently proved 
that below a certain transition temperature, the an- 
ions displace to asymmetric positions, making the two 
TMTTF molecules in a unit cell nonequivalent In 
this Letter, we show that one source of dimerization con- 
trols the existence and the localization length the edge 
states, whereas the other source determines their energy. 

In our theoretical model, we neglect interchain cou- 
pling and consider just one semi-infinite chain. For sim- 
plicity, we will refer to holes as electrons and set h — 1. 
The electron band structure of the chain is described by 
the Hamiltonian 
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where the index n labels the TMTTF sites, s =t, I is the 
spin index, and ipn}s are the creation and annihilation 
operators of electrons. The mean values of the intersite 
tunneling amplitude and the chemical potential are t and 
/i, whereas t and p, are their staggering parts produced by 
dimerization of the chain bonds and sites, respectively. 

The energy dispersion for the Hamiltonian (|^) and (|^) 
is sketched in Fig. ^ For a given chemical potential /z, 
there exist two Fermi points ikp- To study low-energy 
properties of the system, we introduce the operators R 
and L describing the right- and left-moving electrons with 
the momenta close to -\-kp and —kp: 
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FIG. 1. A sketch of the (TMTTF)2X chain. The verti- 
cal rectangles represent the organic molecules TMTTF, and 
the circles indicate the inorganic anions X. The filled objects 
(n — 0, 1, . . .) show the actual chain, whereas the open ob- 
jects (n = 0, —1, —2, . . .) represent an artificial expansion of 
the system by the mirror reflection around the site n = 0. 
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FIG. 2. The energy dispersion for the Hamiltonian 
The single-particle gap at A; = ±7r/2ci opens due to backscat- 
tering ^ caused by the chain dimerization . The correla- 
tion gap at fc = ±7r/4d is caused by umklapp interaction. 
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Here x = nd is the continuous coordinate, and d is the 
average distance between the TMTTF sites. 

The considered system (|^) extends in the positive di- 
rection from the site n = 1 to a length I, which will be 
taken to infinity. The wave functions can be expanded 
to include the site n = with the vanishing boundary 
condition ipQ = 0. Then the system can be artificially 
expanded from the interval (0, 1) to (— Z, I) by the mir- 
ror reflection around the origin n — 0, as shown in Fig. 
|l|, while keeping only the antisymmetric wave functions 
V'n = —ip-n- Because the Hamiltonian couples only 
the nearest neighboring sites, such an expanded prob- 
lem is equivalent to the original semi-infinite problem, 
because there is no direct coupling between the intervals 
{—I, 0) and (0, 1) other than through the point n ~ 0, 
where the wave function vanishes. As Fig. shows, the 
mirror reflection changes the sign of the bond dimeriza- 
tion parameter t, but keeps the same sign for the site 
dimerization p,. Thus, for the expanded problem, Eq. 

should be replaced by t„ — t — sgn(n)(— 1)" f. In 
the continuous expanded system, the boundary condition 
■0(0) = is satisfied by demanding that V'(x) = —ip{—x). 
Using Eq. (0), the latter condition gives 



R{x) = L{-x). 



(4) 



Dimerized chain at half filling. As a warm-up, let us 
study the edge states for noninteracting electrons at half 
filling, where ^ = and kp — n/2d (see Fig. ^). Al- 
though the model and (j^) at half filling does not cor- 
respond to (TMTTF)2X, it does describe some ID poly- 
mers p^ , which have both site and bond dimerization. 
The nondimerized part of the Hamiltonian (|l|) is 



= X! d,^'"F -~Rl{x)idxRs{x) + Ll{x)idxLs{x) 



where vp — 2td is the Fermi velocity. The dimerized part 
of the Hamiltonian gives the backscattering term: 
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The energies e and the wave functions R{x) and L{x) of 
electron eigenstates can be found from the corresponding 
Schrodinger equation 
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Eq. (||) gives the spectrum e = ^Jv'^pk'^ + for the bulk 
states with the energy gap A = 

With the boundary condition '0(0) — 0, which trans- 
lates into -R(O) — L{Q) using Eq. (||), Eq. (^) also gives a 
solution localized at the edge (the Shockley state 



Re{x) = Le{x) 



-t/vF. (7) 



Eq. (0) shows that the energy of the edge state is deter- 
mined by the site dimerization /2, whereas its localization 
length 1/k by the bond dimerization t. The edge state 
exists only when t < 0, i.e. when the edge site (n = 1) 
is connected to the next site (n = 2) by the long bond. 
That configuration is opposite to what is drawn in Fig. 
|l|, where the first bond is short. 

Using Eq. (||), we get Eq. (^) in the expanded picture: 



idxR{x) — [n — itsgn{x)] R{—x) = eR{x). 



(8) 



The edge state (0) maps onto a state localized around 
the kink in the imaginary part of the potential in Eq. 
(^), which is produced by the sign change of the bond 
dimerization t upon mirror reflection around n = 0. Such 
soliton states were studied for polymers in Rcf. . 

Quarter filling. In this case, kp = ir/Ad, as shown in 
Fig. H, and there is no single-particle energy gap at the 
Fermi level. However, a correlation gap (shown as the 
shaded bar in Fig. ||) may develop at the Fermi level due 
to interaction between electrons. Thus, we add the Hub- 
bard interaction U'^^ipl^-^4'ni'4'ni^ni- -^^ ^e rewrit- 
ten in the standard g-o\ogy form in terms of the slow 
fields R and L. Without dimerization, one obtains 

Hint^Y. J dx{2{g^Sss' +gi±Sss')RlLl,R,,L^ 

+2g2RVLl,L^,R^+gi [RIrI,R,,R, + (R ^ L)]} , (9) 

where s =1,1 for s =t, |. 

Introducing the densities of spin and charge 
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the interaction Hamiltonian (|^) can be written as 
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^int = H^:l + h\^\ + 2gi± J2 f RIlIr.L,, (11) 

Ant ^ [ dx [g^pl + pI) - 255 PrPl] , (12) 
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Jo 



(13) 



where ffs = - 2g2- 

r/ie umklapp term. Taking into account that dimer- 
ization (^, treated as a perturbation, hybridizes the 
states with the momenta k = zLkp = zLtt/M and k — 
=F37r/4d, we make the fohowing replacement in Eq. (|^) 



,=pi37™/4 
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where SE = 21^/2. Inserting Eqs. (||) and (|I^ into the 
Hubbard interaction term, we obtain the umklapp term 



^um = Jo dx [g^ + igl) RIr]L .L + H.c 



5,? = -V2Ufid/t, 
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Utd/t. 
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The real part g^ of the umklapp amplitude is propor- 
tional to the site dimerization energy /I, which occurs 
due to the anion displacement, and the imaginary part 
g\ is proportional to the bond dimerization i. 

Bosonization with a boundary. As Eq. (Q) shows, the 
problem on the expanded interval (— /) can be formu- 
lated in terms of the right-moving electrons only. The 
kinetic energy becomes 



Ho — —vf'^^ j dx RI{x) idxR-six) 



(17) 



where vp — td/\/2 for the quarter filling. Imposing the 
periodic boundary condition R{—1) = R{1) together with 
Eq. (^ results in the vanishing boundary condition on 
ip at the other end of the chain: ip{l) = 0. Now we can 
use the standard bosonization technique Q for R. The 
fermion fields Rs (x) can be represented in terms of the 
Bose fields 9s{x): 



R,ix)=f^s e'^=^"^ (20 



-1/2 



(18) 



where C is a short-range cut-off distance of the order of d, 
and fjs are the Majorana fermions: {f7|,?7x} = 0, 77^ = 1. 
The Bose fields obey the commutation relations 

[ds{x),ds' {x')] = in Sss' sgn(x - x'). (19) 

Using the linear combinations 

^P = (^T+^i)/V2, e^^{eT-§^)/V2, (20) 

the charge and spin densities (p^), can be expressed in 
terms of the charge and spin Bose fields (pO|): 



Pr(x) = dx0p{x)/27r, cfb{x) = 5a;(Tp(x)/27r. (21) 

The Hamiltonian ( p^ can be rewritten in terms of the 
densities (M) Rl: 



Ho = nvp I dx [pji{x) + a-j^{x)] 



(22) 



Substituting Eqs. fl) and (^8|) into the interaction Hamil- 
tonian (^ and ( p^ and adding the kinetic energy ([2^), 
we get the full Hamiltonian of the system in the bosonized 
form. It takes the form Hp + H^, where the spin and 
charge sectors are separated. The charge Hamiltonian is 



H„ 



dx [irvp p\{x) - 55 Pr{x)pr{-x) 
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Notice that the imaginary part of the umklapp amplitude 
in Eq. ( p3| ) has a kink, because it is proportional to t ([l^), 
which changes sign between a; < and a; > 0. The spin 
Hamiltonian is 



He = J dx [nva a-%{x) - gi\\aR{x)aR{-x) 
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In Eqs. (H) and (||), Vp^^ =vf± gi/n. 

The quadratic parts of the Hamiltonians (^^ and ( |2^ ) 
can be diagonalized by the Bogoliubov transformations 

()p,cr{x) = Op^a{x) cosh^p^o- - 9p^„{-x) sinh^p^o-, (25) 
Pr{x) = p{x) cosh (f>p + p{~x) sinh (j)p, (26) 
cfr{x) = (7(2;) cosh + (t(— a;) sinh (27) 

where tanh20p = g^/nvp and tanh20o- = gi\\lnva. The 
charge Hamiltonian becomes 



Hr, 



dx 
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iJ2Kp [ep(x)-ep(-x)] 
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(28) 



where Vp = Up/cosh20p and Kp = e^'^f. Similarly, for 
the spin Hamiltonian, we get 



Ha = j dx\Tii„{dxe,{x)/2nf 
(20^ 



(29) 



where = Uo-/ cosh 200- and K^, — e^"^" . 

The umklapp term [the last term in Eq. (p8|)] is irrele- 
vant when Kp > 1 H. In the opposite case Kp < 1, this 



3 



term is relevant and leads to opening of an energy gap 
in the charge sector. At the particular value Kp = 1/2, 
the charge Hamiltonian Hp can be mapped to a model 
of noninteracting fermions with an energy gap ||7|,|l3|| . By 
introducing the effective fermions (holons) 



V'p(a;) 



(2^0 
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and taking into account the commutation relations (19) 
the last term in Eq. (28) can be written as 



il0p{x)-epi-x)] 



-i sign(x) 



Then the charge Hamiltonian (Esh becomes 



Hp= I dx ^^-Vp^pl{x)i^x'^ijp{a 



(31) 



[53' - z<?3« sgn(x)] ^Plix)^Pp{-x)/27rC} . (32) 



The energy spectrum Sp of the Hamiltonian ( |32| ) can 
found from the corresponding Schrodinger equation 



ivpdx'ippix) 



93 - sgn(a;) - 



27rC 



i^pi-x) ^ EpTppix). (33) 



Eq. (^) is formally equivalent to Eq. (^ with and 
playing the roles of /2 and i, respectively. However, as 
Eq. dig ) shows, (x t and 173 oc p,, so the roles of jl and 
t are reversed in Eq. ( ^3|) compared with Eq. ( p|). This 
is a consequence of the factor — isign(x) in Eq. (|3l|). Eq. 

( ^ ) gives the spectrum Sp — yjvpk'^ + Ap for the bulk 

states with the charge gap Ap = ^/{gf)^+{g^ /2Tr(. It 
also gives an edge solution localized around x = 0: 



-pe 



-gy27rC = -Apt/^P + 2fl^, 



(34) 
(35) 
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where we used Eq. (|6|) assuming f7 > for the repulsive 
Hubbard model. We see that, opposite to the electron 
edge state (^, the energy (^5|) of the holon edge state 
( ^ ) is determined by the bond dimerization t and its 
localization length ( ^6|) is controlled by the site dimer- 
ization p.. The holon edge state exists only when 173 < 0, 
i.e. p, > 0. Taking into account Eq. (||), we conclude 
that the holon edge state exists when the chemical po- 
tential of the edge site {n = 1) is higher, i.e. when the 
anion displaces toward the edge site, as shown in Fig. 
|l|. Thus, we predict that the holon edges states should 
appear in (TMTTF)2X below the anion transition tem- 
perature g. Their energies (pq) depend on whether the 
first bond at the edge is short or long. 

Similar consideration S can be done for the spin 



Kcr < 1, i.e. < 0, whereas the required condition 
for the existence of an edge state is gi± > 0. However, 
gi± = 9i\\ ill a physical system because of the spin SU(2) 
invariance. So, it is not possible to satisfy the conditions 
for opening of a gap and existence of an edge state simul- 
taneously. Thus, we conclude, contrary to Ref. that 
spinon edge states cannot exist. 

In conclusion, we have shown that holon edge states 
may appear in a ID model of interacting electrons with 
an energy gap in the charge sector. The existence of 
the holon edge state and its localization length are con- 
trolled by the real part of the umklapp interaction am- 
plitude, whereas its energy is determined by the imagi- 
nary part. In the quasi-one-dimensional organic conduc- 
tors (TMTTF)2X, these umklapp amplitudes are propor- 
tional to the site and bond dimerizations, correspond- 
ingly. The bond dimerization is always present in these 
materials. However, the site dimerization occurs only be- 
low the recently discovered ferroelectric transition, where 
the inorganic anions X displace to asymmetric positions 
|]lT| . We predict that the holon edge states should ap- 
pear below the anion transition. They could be probed 
by electron tunneling into the ends of ID chains or by 
other types of spectroscopy. 
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